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BACKGROUND  AND  OBJECTIVES 


In  July  1976,  a  contract  was  awarded! by  the  Office  of  Naval 
Research  to  Electric  Boat  Division  iri  the  area  of  computer 
analysis  of  welded  ship  structures .*<7 The  overall  objective 
of  this  research  was  to  apply  and  further  develop  where 
necessary  the  technology  for  predicting  distortion  of  sub¬ 
marine  structure  due  to  welding.  The  underlying  motivation 
in  providing  this  technology  is  that  it  can  be  used  in  design 
and  construction  of  submarines  to  help  determine  the  most 
cost  effective  weld  joint  designs  and  welding  procedures. 

The  specific  objectives  of  the  initial  contract  were: 

/ 

4)  implementation  of  MIT  finite  element  computer 

programs  for  prediction  of  temperatures,  residual 
stresses,  and  distortion  due  to  welding^ 

r 

2)  evaluation  of  MIT  computer  programs  on  a  fusion 
weld  of  two  plates  and  comparison  with  results 
obtained  by  MARC-CDC  finite  element  program^ 

3)  analysis  of  a  butt  welded  unrestrained  plate 
with  MIT  computer  codes  and  comparison  with 


test  results. 


The  results  of  that  study  are  documented  in  Reference  1 . 

In  June  1977,  a  contract  modification  extended  the  scope  of 
the  initial  contract  to  the  following  major  technical 
objectives : 

1)  transient  thermoplastic  analysis  of  a  butt 
welded  unrestrained  plate  using  the  Electric 
Boat  Division  modified  NONSAP  finite  element 
computer  code, 

2)  further  development  and  application  of  the 
Shrinkage  Force  Method  (SFM)  to  welding 
analysis  including  investigation  of  weld 
bead  sequencing,  non-symmetry  in  the  weld 
groove  and  remelt  effects.  The  SFM  is  a 
simplified  analytical  approach  for 
predicting  weld  distortion. 

The  results  of  that  study  are  documented  in  Reference  2. 
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A  second  contract  modification  was  awarded  in  September  1978 
to  apply  the  developed  analytical  methods  to  circumferen¬ 
tial  butt  welding.  An  experiment  conducted  at  the  Massachu¬ 
setts  Institute  of  Technology  on  an  HY-130  cylinder  butt 
weld  (Reference  3)  was  chosen  as  the  test  problem  for  further 
evaluation  of  the  transient  thermoplastic  analysis  and 
Shrinkage  Force  Method  developed  under  the  initial  contract 
and  first  modification.  In  addition  application  of  the 
Shrinkage  Force  Method  to  multipass,  circumferential  welding 
of  submarine  HY-80  ring  stiffened  cylinder  sections  was 
also  included  in  the  workscope.  The  analytical  approach 
employed  in  this  latter  task  was  developed  under  Electric 
Boat  Division  IRAD  Task  573.  It  utilizes  a  layered  shell 
approximation  to  simulate  the  weld  build-up. 

The  results  of  these  latest  studies  are  documented  in  this 
report . 


II.  THERMOMECHANICAL  DISTORTION  ANALYSIS  OF  MULTIPASS 

GIRTH  WELDING  IN  AN  HY-130  CYLINDER 

The  objective  of  this  segment  of  the  contract  was  to  analy¬ 
tically  predict  weld  distortion  due  to  girth  welding  of  an 
HY-130  cylinder  and  compare  the  results  with  available  test 
data.  The  cylinder  considered  was  the  18"  diameter,  HY-130 
test  sample  (see  Figure  II-C-2)  GMAW  welded  at  the  Massa¬ 
chusetts  Institute  of  Technology  and  reported  in  Reference  3. 

The  approach  to  the  mathematical  modelling  problem  was  as 
follows  s 

a.  Analytically  predict  the  transient  temperature 
distribution  during  girth  welding. 

b.  Analytically  predict  transient  thermoplastic 
distortions  using  the  above  predicted  transient 
temperatures. 

c.  In  addition  apply  the  Shrinkage  Force  Method 
to  predict  thermomechanical  distortions. 

Detailed  discussions  of  the  analyses  and  results  are  presented 
in  the  following  sections. 
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II-A  HEAT  TRANSFER  ANALYSIS 


o  Introduction 

A  heat  transfer  analysis  was  made  using  the  TEMP  computer 
program,  in  order  to  calculate  the  temperature  transients 
required  by  the  ensuing  NONSAP  structural  analysis.  This 
thermal  transient  program,  which  uses  the  finite  element  analy¬ 
sis  technique,  is  based  on  the  work  described  in  Reference  4. 
The  program  was  obtained  in  1978  during  a  computer  workshop 
at  Union  College  in  Schenectady,  New  York. 

The  program  can  be  used  on  either  one  or  two  dimensional  plane 
or  axisymmetric  bodies.  Boundary  conditions  of  convection, 
radiation,  heat  flux  or  nodal  temperature  may  be  specified, 
and  all  boundary  conditions  can  be  functions  of  time. 

Five  different  materials  ma”  be  used  with  the  associated 
thermal  properties  considered  to  be  quadratic  functions  of 
temperature  about  a  specified  reference  temperature.  The 
proper  selection  of  coefficients  and  reference  temperature 
can  yield  either  a  linear  property  variation  with  temperature 
or  a  constant  property  value. 
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The  utilization  of  the  TEMP  computer  program  also  allowed  the 
use  of  the  same  model  elemental  break-up  as  used  by  NONSAP. 
Thus,  the  required  nodal  temperatures  were  computed  directly, 
and  the  need  for  spatial  interpolation  was  eliminated. 


The  analytical  model,  excluding  the  elemental  break-up,  is 
shown  in  Figure  II-A-1.  The  model  is  a  two  dimensional,  axi- 
symmetric  representation  of  the  HY-130  cylinder.  In  this  case 
the  TEMP  program  uses  elemental  depths  (lengths  in  the  circum¬ 
ferential  direction)  associated  with  an  angle  of  one  radian. 

This  yields  element  depths  of  approximately  nine  inches.  The 
weld  region  is  represented  as  three  weld  layers  rather  than 
the  actual  six  bead  lay-up  which  is  somewhat  questionable  as 
to  its  thermal  modelling  accuracy,  since- each  layer  is  composed 
of  several  beads  which  are  actually  deposited  at  different  times. 
(See  Figure  II-C-10.)  The  elements  and  nodes  associated  with 
the  three  weld  layers  are  shown  in  Figure  II-B-2.  The  first 
weld  layer  is  composed  of  elements  6001  through  6009,  the 
second  weld  layer  is  composed  of  elements  6010  through  6015, 
and  the  third  weld  layer  is  composed  of  elements  6016  through 
6021.  Each  weld  layer  deposition  was  modelled  as  an  instan¬ 
taneous  addition  of  the  weld  metal  at  its  vaporization  temperature 
of  5000°F  at  time  equal  to  zero  resulting  in  a  separate  thermal 


analysis  for  each  of  the  three  weld  layer  depositions.  It 
was  also  assumed  that  the  base  metal  and  any  previously  deposited 
weld  layers  were  preheated  to  a  temperature  of  200°F  at  time 
equal  to  zero.  The  resulting  initial  condition  for  each  tem¬ 
perature  transient  calculation  was  that  all  nodes  of  the  weld 
layer  being  deposited  were  at  a  temperature  of  5000°F,  and 
all  other  nodes  of  the  structure  were  at  a  temperature  of 
200°F .  As  the  transient  progressed  the  deposited  weld  layer, 
which  was  assumed  to  have  the  same  thermal  properties  as  HY-130, 
was  cooled  by  conduction  to  the  rest  of  the  structure  since 
it  was  assumed  that  the  convective  and  radiative  heat  losses 
were  negligible  as  indicated  by  the  insulated  surfaces  on 
Figure  II-A-1. 

The  heat  addition  approach  used  above  differs  from  the  "ramp 
heat  flux"  heat  addition  approach  previously  used  for  flat 
plate  welding  thermal  analysis  (Reference  1) .  The  ramp  approach 
assumes  an  applied  heat  flux  which  starts  at  a  value  of  zero 
at  zero  time,  increases  to  a  prescribed  value  at  a  specified 
time,  continues  at  a  constant  value  for  a  given  period  of 
time,  then  decreases  to  zero  at  a  certain  time.  It  was  felt 
that  the  ramp  heat  addition  approach  would  not  offer  a  good 


representation  of  the  actual  heat  input  for  this  case,  since 
the  ramp  would  have  to  extend  over  the  first  45  seconds  of 
the  thermal  transient. 

It  was  determined  that  the  total  heat  released  by  the  weld 
metal  during  a  cooldown  from  5000°F  to  200°F  would  be  939  Btu 
per  pound.  Since  the  latent  heat  of  fusion  of  117  Btu  per 
pound  would  be  a  small  portion  of  the  total  (approximately 
12%),  it  was  neglected  in  order  to  simplify  the  analysis. 

The  thermal  properties  for  HY-130  steel  are  presented  in 
Table  II— A- 1 .  The  following  method  was  used  to  account  for 
the  thermal  property  variations,  since  the  variations  do  not 
lend  themselves  to  either  the  quadratic  or  linear  represen¬ 
tations  of  the  TEMP  computer  program.  The  structure  was 
divided  into  two  sections;  namely,  a  large  low  temperature 
(less  than  or  equal  to  400°F)  section,  and  a  small  high 
temperature  (greater  than  400°F)  section.  The  two  sections 
were  composed  of  different  materials  with  different  constant 
values  for  their  thermal  properties.  The  properties  for  each 
section  were  evaluated  by  taking  an  integrated  average  over 
the  appropriate  temperature  range. 
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o  Results 


Typical  predicted  temperature  transients  for  the  first,  second 
and  third  weld  layer  depositions  for  four  nodal  points  are 
presented  in  Figures  II-A-2  through  II-A-4.  The  nodal  point 
locations  were  selected  such  that  the  predictions  for  the 
first  weld  layer  deposition  could  be  compared  to  the  experi¬ 
mental  data  of  Reference  3.  in  Reference  3,  experimental  data 
is  presented  for  the  first  weld  pass  only.  It  is  unfortunate 
that  experimental  temperature  data  is  not  available  for  the 
other  weld  passes  also.  These  comparisons  are  presented  in 
Figures  II-A-5  through  II-A-8.  As  indicated  by  the  figures, 
the  predicted  temperatures  for  the  first  weld  layer  deposition 
are  substantially  lower  than  the  experimental  values.  It  is 
felt  that  this  discrepancy  could  be  due  to  the  use  of  inaccurate 
low  values  of  specific  heat  capacity  for  the  liquid  metal  weld 
layer  at  early  times  of  the  predicted  transient.  A  substan¬ 
tially  higher  value  would  yield  a  higher  temperature  differ¬ 
ence  driving  force  for  conduction,  as  well  as  more  energy 
available  for  base  metal  heat-up. 
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THERMAL  ANALYTICAL  MODEL 


FIGURE 


PREDICTED  TEMPERATURE  TRANSIENTS  FOR 
POINTS  LOCATED  ON  OUTER  SURFACE  OF 
CYLINDER  FOR  FIRST  WELD  LAYER  DEPOSITION 


<j0  (N.  FROM  WELD  CENTERLINE 
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FIGURE  I- A- 2 


SECONDS 


FIGURE  H-A-3 


PREDICTED  TEMPERATURE  TRANSIENTS  FOR 
POINTS  LOCATED  ON  OUTER  SURFACE  OF 
CYLINDER  FOR  THIRD  WELD  LAYER  DEPOSITION 


1.0  IN.  FROM  WELD  CENTERLINE 


TEMPERATURE 


COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED 
TEMPERATURE  TRANSIENT  AT  A  POINT 
LOCATED  1.5  IN.  FROM  WELD  CENTERLINE 
ON  OUTER  SURFACE  OF  CYLINDER 
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FIGURE  H-A-6 


COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED 
TEMPERATURE  TRANSIENT  AT  A  POINT 
LOCATED  2.0  IN.  FROM  WELD  CENTERLINE 
ON  OUTER  SURFACE  OF  CYLINDER 
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FIGURE  H-A-7 


TEMPERATURE 


comparison  of  experimental  and  predicted 

TEMPERATURE  TRANSIENT  AT  A  POINT 
LOCATED  2.5  IN.  FROM  WELD  CENTERLINE 
ON  OUTER  SURFACE  OF  CYLINDER 
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FIGURE  H-A-8 


II-B  THERMOPLASTIC  ANALYSIS  (USING  THE  NONSAP 
COMPUTER  PROGRAM) 

o  Introduction 

In  1977  the  NONSAP  finite  element  program  was  modified  to 
perform  thermoplastic  analysis  with  temperature  dependent 
material  properties.  Additional  modifications  were  made  in 
1978  to  assist  in  the  analysis  of  weld  distortion  problems. 
Specifically,  a  "birth"  and  "death"  option  was  added  for 
two  dimensional  finite  elements  that  represent  the  weld 
deposit.  At  specified  times  during  the  analysis,  the  ele¬ 
ments  in  the  weld  groove  region  may  be  activated  (achieve 
strength  from  cooldown)  or  deactivated  (melt) .  This  added 
feature  to  the  program  greatly  assists  in  the  modelling  of 
the  weld  region. 

In  1979  the  modified  NONSAP  program  was  applied  to  the  HY-130 
cylinder  girth  weld  reported  in  Reference  3.  The  transient 
temperature  distributions  used  in  this  analysis  were  cal¬ 
culated  via  the  TEMP  computer  code  as  reported  above  in 
Section  II-A. 
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o  Analysis  Methods 


Figure  II-B-1  shows  the  overall  NONSAP  finite  element  model 
used  for  the  thermoplastic  analysis.  Figure  II-B-2  shows 
the  model  of  the  weld  groove  region  which  is  a  three  layer 
representation  of  the  six  weld  beads  used  in  the  actual 
cylinder  welding.  Figures  II-B-3  thru  II-B-7  show  the  re¬ 
mainder  of  the  model  elements.  The  model  employed  symmetric 
boundary  conditions  in  the  weld  groove  region  and  free  move¬ 
ment  at  the  opposite  end  of  the  cylinder.  The  "+y"  direction 
of  the  local  coordinate  system  was  in  the  radial  direction 
and  "+z"  ran  along  the  length  of  the  shell. 

Four-noded  isoparametric,  quadrilateral,  plain-strain  finite 
elements,  with  a  2x2  Gaussian  integration  scheme,  were  util¬ 
ized  in  the  analysis.  The  elements  were  reduced  to  triangular 
shape  where  necessary,  to  represent  the  geometry  and  to  make 
transitions  from  large  element  size  to  small  element  size 
as  the  weld  region  is  approached. 

From  the  transient  heat  transfer  analysis,  a  file  of  nodal 
temperature  distributions  as  a  function  of  time  were  generated 
to  simulate  the  heatup-cooldown  cycle  caused  by  individual 
deposition  of  the  three  weld  layers.  The  temperature  dependency 
of  the  material  properties  used  in  this  analysis  are  plotted 
in  Figure  II-C-1. 
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The  math  model  used  for  the  heat  transfer  analysis  was 
identical  to  the  NONSAP  analysis  math  model.  This  allowed 
a  sequence  of  temperature  distributions  from  time  =  0.0 
seconds  to  time  =  60.0  seconds.  The  maximum  temperature 
change  between  any  two  nodes  in  going  from  time  (i)  to  time 
(i  +  1)  ranged  from  0.0  to  4800°F. 

From  these  nodal  temperatures,  the  temperature  at  the  four 
Gauss  points  in  each  finite  element  is  calculated  in  NON SAP . 

The  tangent  stiffness  is  then  calculated  using  the  material 
properties  associated  with  the  current  temperature  and  stresses 
at  the  Gauss  points.  The  increment  of  thermal  load  in  going 
from  time  (i)  to  time  (i  +  1)  is  applied  to  obtain  the  incre¬ 
ment  in  nodal  deflections  and  the  associated  incremental 
stresses  at  the  Gauss  points  assuming  elastic  behavior. 

Then,  the  total  stress  at  each  Gauss  point  is  compared  to  the 
Von  Mises  yield  surface  at  the  current  temperature  to  determine 
whether  the  increment  in  stress  calculated  elastically  is 
admissable  or  whether  the  yield  surface  has  been  exceeded. 

If  the  yield  surface  has  been  exceeded,  the  stresses  are 
adjusted  to  be  consistent  with  the  specified  plastic  behavior 
model. 

Equilibrium  iterations  in  the  NONSAP  program  are  performed 
until  the  displacement  has  converged  to  within  a  0.005  error 
tolerance.  Within  each  of  these  iterations,  the  stresses 
are  updated,  checked  against  the  yield  surface,  and  corrected 
if  necessary. 
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o  Analysis  Results 

The  initial  approach  was  to  analyze  each  of  the  weld  layers 
vrith  a  successive  series  of  heatup-cooldown  cycles  intro¬ 
ducing  weld  layers  1  thru  3  at  the  appropriate  times.  The 
approach  had  to  be  abandoned  during  the  cooldown  of  weld 
layer  1  (elements  6001  to  6009)  ,  when  unresolvable  conver¬ 
gence  difficulties  developed.  A  considerable  effort  was 
made  to  obtain  convergence  of  the  solution,  however  no 
solution  strategy  was  found  to  accomplish  this.  When  the 
solution  was  continued  to  the  next  time  step,  the  displace¬ 
ment  solution  blew  up.  The  following  is  a  possible  reason 
the  analysis  would  not  converge: 

The  nodal  temperature  at  the  four  Gauss  points  in  each  finite 
element  is  calculated.  The  tangent  stiffness  of  the  element 
is  then  calculated  using  the  material  properties  associated 
with  the  current  temperature.  The  weld  region  which  has 
small  elements  experiences  large  stiffness  changes  due  to 
the  temperature  transients  of  the  heatup  and  cooldown  cycles. 
The  remaining  portion  of  the  model  which  has  relatively 
larger  elements  maintains  an  average  stiffness  due  to  the 
smaller  changes  in  temperature.  The  combination  of  the  weld 
region  experiencing  drastic  changes  in  stiffness  and  the 
larger  elements  having  a  relatively  constant  stiffness  may 
be  causing  an  ill-conditioned  matrix. 
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THERMOPLASTIC  MATH  MODEL 

FIGURE  II- B- I 
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FIGURE  H-  B  -  2 
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II-C.  SHRINKAGE  FORCE  METHOD 


o  Introduction 

This  weld  distortion  prediction  technique  is  based  upon  the 
assumption  that  the  shrinkage  forces  developed  in  each  weld 
layer  during  thermal  cooldown  are  the  major  contributors  to 
distortion.  In  the  application  of  the  method,  the  geometry 
of  each  weld  layer  along  with  the  appropriate  temperature 
changes  (AT)  for  the  cooldown  of  the  weld  metal  must  be  known. 
In  its  simplest  form  the  linear  coefficient  of  thermal 
expansion  (aT)  and  Young's  modulus  (E)  are  used  to  calculate 
the  linear  contraction  force  layer  by  layer.  This  is  applied 
to  the  workpiece  in  a  series  of  numerical  solutions  for  dis¬ 
tortions  .  The  procedure  can  approximately  account  for  the 
variation  with  temperature  of  material  properties  (E  and  a 
in  the  elastic  case)  by  simply  using  an  integrated  average 
over  the  cooldown  temperature  range.  The  primary  feature 
of  the  SFM  method  is  the  elimination  of  the  need  for  an 
expensive  time  dependent  nonlinear  heat  transfer  analysis. 
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o  Idealization  of  Material  Properties 


The  material  properties  used  in  applying  the  Shrinkage  Force 
Method  of  analysis  to  the  NONSAP  computer  code  are  presented 
in  Table  II-C-1  and  Figure  II-C-1.  NONSAP  is  limited  to 
six  (6)  data  points  to  define  material  properties  as  a  func¬ 
tion  of  temperature  and  assumes  a  linear  relationship  between 
the  data  points.  A  temperature  of  2500°F  was  used  as  the 
starting  point  for  cooldown  in  the  analysis  since  this  point 
represents  the  temperature  where  the  weld  metal  starts  to 
achieve  strength. 

o  Finite  Element  Model 

The  finite  element  model  developed  and  used  to  predict  dis¬ 
tortion  caused  by  multipass  welding  in  the  HY-130  cylinder 
is  shown  in  Figure  II-C-3.  A  twenty  inch  long  segment  of 
the  eighteen  inch  diameter  cylinder  was  modeled  with  a  60 
degree  included  angle  single  vee  outside  weld  prep  with 
1/8  inch  land.  The  weld  region  consisted  of  a  three  layer 
representation  of  the  six  weld  beads  used  to  fill  the  joint 
shown  in  Figure  II-C-2.  A  detailed  breakup  of  the  finite 
element  model  showing  joint  and  element  numbering  is  shown 
in  Figures  II-C-4  thru  II-C-10.  This  model  is  a  mirror 
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image  of  the  one  used  in  the  previous  section.  The  origin 
of  the  model  was  located  at  the  weld  centerline  with  posi¬ 
tive  "y"  along  the  radial  direction  and  position  "z"  along 
the  length  of  the  cylinder.  Symmetry  boundary  conditions 
were  applied  at  the  weld  centerline  (0,0)  while  the  opposite 
edge  of  the  cylinder  was  free. 

o  Analytical  Approach 

Similar  to  the  approach  used  in  the  previous  section,  the 
SFM  idealized  the  cylinder  as  a  two  dimensional  axisymmetric 
problem.  The  model  employed  two  dimensional  four  noded  iso¬ 
parametric,  quadrilateral,  plain  strain  finite  elements, 
with  a  2x2  Gaussian  integration  scheme.  The  elements  were 
reduced  to  triangular  shapes  where  necessary  to  represent 
the  geometry  and  for  transition  from  large  element  size  to 
small  element  size  near  the  weld  region. 

The  HY-130  material  was  assumed  to  have  six  different  moduli 
(E  =  1,2, 3,... 6)  and  thermal  coefficients  (aT  =  1,2, 3,... 6) 
in  each  layer  analysis.  Values  for  E,  aT,  and  AT  are  given 
in  Table  II-C-1  and  shown  in  Figure  II-C-1. 
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The  distortion  for  each  weld  layer  as  it  cools  down  is 
computed  and  displacements  are  simply  accumulated  layer  by 
layer.  The  model  reflects  the  changing  geometry  as  each 
layer  is  added.  Internal  stresses  are  also  accumulated 
throughout  the  analysis.  When  the  weld  elements  are  acti¬ 
vated  at  T  =  2500°F,  the  internal  stresses  within  the  pre¬ 
viously  cooled  material,  which  is  then  assumed  to  remelt, 
are  set  equal  to  zero  (i.e.  no  strength  is  assumed),  and 
the  change  in  deflection  is  computed  for  this  condition. 

The  deflections  and  internal  stresses  are  accumulated  for 
solidification  and  remelt  cycles  until  the  final  deformed 
state  is  determined. 

Equilibrium  iterations  are  performed  until  the  displacement 
solution  has  converged  to  within  a  0.005  error  tolerance. 

Within  each  iteration,  the  stresses  are  updated. 

The  approach  used  to  analyze  the  HY-130  cylinders  was  as 
follows : 

1)  Start  heatup  cycle  for  layer  1  at  200°F  (preheat  TIME=0.0 

temperature)  for  all  joints  in  the  model. 

2)  Heat  joints  along  the  weld  prep's  edge  to 
2500°F,  to  make  them  compatible  with  the 
layer  1  weld  joints  to  be  introduced. 


TIME=1 . 0 
to 

TIME=10 . 0 
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3)  Activate  layer  1  weld  joints  and  elements 


TIME  =11.0 


of  the  model  at  a  temperature  of  2500°F. 

TIME  =12.0 
to 

TIME  =20.0 

5)  Begin  heatup  cycle  for  layer  2.  A  problem  was  TIME  =21.0 

encountered  at  this  point  in  re-starting  the 
analysis  for  second  layer  heatup;  (T  =  500°F) 
due  to  no  convergence  at  affected  nodes. 

The  solution  to  this  problem  was  to  go  up  in  small  temperature 
increments  at  first,  until  good  convergence  was  reached.  The 
following  incremental  runs  were  used  to  reach  500°F. 

TIME _ TEMP.  AT  AFFECTED  NODES  °F 

20.2  225 

20.4  250 

20.6  300 

20.8  400 

21.0  500 

Convergence  was  reached  at  time  =  21.0,  T  =  500°F  and  saved. 


4)  Cooldown  layer  1  model  from  2250°F  to  200°F, 
preheat  temperature. 


Continued  heatup  for  layer  2. 


TIME  =22.0 
TO 

TIME  =29.0 


6) 


TIME=30.0 


Activate  layer  2  weld  elements  and  joints  at 
T  =  2500°F . 


7)  Start  cooldown  cycle  for  layer  2  from  2250°F  to 
200°F.  Again  a  problem  was  encountered  due  to  no 
convergence-  The  solution  to  this  was  to  let  the 
temperature  stay  at  2500°F  for  two  analyses  to 
achieve  good  convergence;  then  start  dropping 
temperature  with  small  AT's  to  2250°F.  The  fol¬ 
lowing  runs  were  made ; 


TIME  TEMP.  AT  AFFECTED  NODES  °F 


30.2 

2500.0 

30.4 

2500.0 

30.6 

2475.0 

30.8 

2400.0 

31.0 

2250.0 

Good  convergence  was  achieved  and  the  analysis  proceeded 
through  the  cooldown  cycle. 

Continue  cooldown  to  preheat  temperature  for  all 
affected  nodes. 


TIME=31 . 0 


TIME=32 . 0 
TO 


TIME=40 . 0 


NOTE 


Due  to  the  convergence  difficulties  encountered 
in  the  heatup  and  cooldown  cycles  for  layer  two, 
interim  analyses  were  made  for  layer  three  to 
avoid  these  problems  again. 

8)  Heatup  for  layer  3. 

9)  Activate  layer  3  weld  joints  and  elements 
at  T  =  2500°F . 

10)  Cooldown  cycle  for  layer  3. 

11)  Analysis  complete. 
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TIME=4 1 . 0 
TO 

TIME=49 . 0 


TIME=50 . 0 


TIME=51 . 0 
TO 

TIME=60 . 0 


V  S 


o  Analysis  Results 


Figure  II-C-11  shows  the  cumulative  radial  distortion  at  the 
weld  centerline  versus  the  NONSAP  run  sequence  number.  The 
run  sequence  number  corresponds  to  an  applied  material  temper¬ 
ature  (°F)  used  in  the  NONSAP  analysis.  These  are  tabulated 
in  Table  II-C-2  along  with  the  corresponding  temperature  and 
heatup  or  cooldown  cycles  fnr  each  layer.  The  distortion 
pattern  for  each  layer  is  generally  characterized  by  distortion 
relief  during  the  heatup  cycle  and  increasing  distortion  during 
the  cooldown  cycle. 

Figure  II-C-12  shows  the  layer  by  layer  radial  distortion 
along  the  length  of  the  cylinder.  As  can  be  seen  the  greatest 
amount  of  distortion  occurs  within  2.0  inches  from  the  weld 
centerline  for  each  layer,  and  rapidly  diminishes  after  that 
point.  A  final  distortion  prediction  of  -0.0194  inches  is 
shown  at  the  weld  centerline. 


j 
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RUN 

TIME" 

TEMP 

Op 

THERMAL  CYCLE 

0.0 

200° 

1.0 

500° 

2.0 

800° 

3.0 

1000° 

4.0 

1200° 

5.0 

1400° 

6.0 

1600° 

7.0 

1800° 

8.0 

2000° 

9.0 

2250° 

10.0 

2500°F 

11.0 

2250° 

12.0 

2000 

13.0 

1800° 

14.0 

1600° 

15.0 

1400 

16.0 

1200° 

17.0 

1000° 

18.0 

800° 

19.0 

500° 

20.0 

200 

21.0 

500° 

22.0 

800° 

23.0 

1000° 

24.0 

1200 

25.0 

1400° 

26.0 

1600° 

27.0 

1800° 

28.0 

2000° 

29.0 

2250° 

30.0 

2500° 

31.0 

2250° 

32.0 

2000° 

33.0 

1800° 

34.0 

1600° 

35.0 

1400° 

36.0 

1200° 

37.0 

1000 

38.0 

800° 

39.0 

500° 

40.0 

200° 

Layer  1  heatup 


End  Activate  Layer 

1 


End 

Layer  2  heatup 


End  Activate  Layer 

Layer  2  cooldown  2 


End 


TABLE  II-C-2 
NONSAP  RUN  "TIMES" 
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41.0 

500° 

42.0 

80g° 

43.0 

100° 

44.0 

1200° 

45.0 

1400° 

46.0 

1600° 

47.0 

1800° 

48.0 

2000° 

49.0 

2250° 

50.0 

2500° 

51.0 

2250° 

52.0 

2000° 

53.0 

1800° 

54.0 

1600° 

55.0 

1400° 

56.0 

1200° 

57.0 

1000° 

58.0 

800° 

59.0 

500° 

60.0 

200° 

Layer  ^  heatup 


End  Activate  Layer 

Layer  3  cooldown  3 


End 


End  of  Analysis 


TABLE  II-C-2 
(CONTINUED) 
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WELD  JOINT  GEOMETRY  FOR  HY*  1 30  CYLINDER 
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FIGURE  E-C-2 
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PLANE  OF  SYMMETRY - 1 


NONSAP  MATH  MODEL 
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FIGURE  H-C-3 


aw. 


1 

j 


?1 

?a 

73 

74 

75 

76 

77 

78 

79 


NONSAP  MATH  MODEL 
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FIGURE  H-C-7 
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FIGURE  H-C-8 


_ A-ii? 


^  S029 

5028 

—i 

U±> 

_^_L*L - 

rn 


-  A 

--  A 


138 

131 

132 


NONSAP  MATH  MODEL 


49 


FIGURE  II- C- 9 


RADIAL  DEFLECTION  (*  10  INCHES) 


NOTE:  AFTER  LAYER  I  COOLDOWN  ’0.0054 

AFTER  LAYER  2  COOLDOWN  ’0.0129 

AFTER  LAYER  3  COOLDOWN  ’0.0194 


□  LAYER  I 
O  LAYER  2 
A  LAYER  3 


WELD 

L 


INCREMENTAL  LAYER  RADIAL  DEFLECTIONS 
vs.  CYLINDER  LENGTH 
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FIGURE  E-C-12 


III.  SHRINKAGE  FORCE  METHOD  ANALYSIS  OF  SUBMARINE 
HULL  CYLINDER  GI RTH  WELD 

o  Background 

In  1977,  the  hull  butt  between  cylinders  20  and  21  of  the 
SSN704  was  selected  by  the  Welding  Engineering  Department 
for  evaluation  of  a  single  vee  hull  butt  joint  design,  with 
ceramic  backing  strip  on  the  outside  of  the  hull,  and  all 
welding  performed  from  the  inside  of  the  hull.  The  primary 
objective  of  the  Welding  Engineering  test  was  to  see  if  it 
was  possible  to  minimize,  if  not  completely  eliminate  back- 
gouging.  Special  steps  were  taken  to  control  the  weld  prep 
geometry  and  to  carefully  monitor  the  weld  distortion.  Since 
the  weld  prep  geometry  and  actual  welding  of  the  joint  was 
tightly  controlled,  it  provided  an  excellent  opportunity  to 
obtain  extensive  distortion  measurements.  These  experimental 
measurements  have  proven  to  be  an  excellent  data  base  for 
evaluation  of  analytical  methods  for  predicting  weld  distor¬ 
tion  of  hull  butt  joints.  The  detailed  distortion  measurements 
and  a  description  of  the  measurement  procedure  are  given  in 
Reference  6. 
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The  original  analytical  approach,  in  which  the  molten  metal 
was  assumed  to  behave  elastically,  was  applied  to  this  pro¬ 
blem  in  1977,  using  the  B0S0R4  computer  code  (see  Refer¬ 
ence  7) .  The  conservativeness  of  this  approach  led  to 
an  overprediction  of  final  distortion  by  a  factor  of  approx¬ 
imately  2.  In  1978,  this  joint  was  re-analyzed  using  the 
B0S0R5  code  (see  Reference  8)  to  include  plastic  material 
behavior  of  the  weld  metal.  The  main  difference  in  this 
approach  is  that  the  weld  metal  during  cooldown  is  limited 
to  developing  a  maximum  Von-Mises  effective  stress  equal 
to  the  specified  filler  metal  yield  strength  at  room 
temperature.  This  change  provides  an  upper  limit  on  the 
maximum  force  that  the  cooling,  shrinking  weld  metal  can 
exert  on  the  pieces  being  welded.  The  distortion  of  the 
elastic-plastic  approach  was  approximately  50  percent  of 
the  previous  analysis  and  quite  closely  approximates  the 
experimental  results.  See  Figure  III-l. 

Improvement  of  the  analytical  approach,  the  "shrinkage 
force"  method,  was  undertaken  in  1979.  The  effort  concen¬ 
trated  on  two  areas:  one  was  to  gather  additional  experi¬ 
mental  data  on  the  amount  of  remelt  for  various  welding 
techniques.  The  other  was  to  include  the  effects  of  remelt 
in  the  analytical  model.  The  results  of  this  effort  are 
reported  in  this  section. 


o  Experimental  Results 

The  extent  of  base  plate  melt  which  occurs  during  welding 
had  previously  been  assumed  to  be  1/8"  on  each  side  of  the 
weld  joint.  This  assumption  was  based  on  the  judgement  of 
welding  engineers.  This  melted  base  plate  was  assumed  to 
act  with  the  deposited  weld  metal,  thus  increasing  the  trans¬ 
verse  distance  across  the  weld  joint  by  1/4"  in  the  analytical 
model  of  the  weld  region.  Since  weld  test  samples  for  a 
number  of  different  weld  processes  are  now  available,  it 
was  possible  to  evaluate  this  assumption  quantitatively. 

The  results  showed  that  base  plate  melt  is  significantly 
less  than  originally  assumed  -  ranging  from  1/16"  to  1/64", 
and  generally  about  1/32"  (see  Table  III-l) .  This  is  attributed 
to  the  fact  that  the  plates  being  welded  provide  a  large 
heat  sink  to  rapidly  conduct  the  welding  heat  away  from  the 
sides  of  the  weld  joint,  thus  minimizing  the  lateral  base 
plate  melt. 

The  experimental  data  gathered  also  shows  that  the  only 
significant  melting  which  occurs  as  a  result  of  depositing 
a  weld  bead  is  remelt  of  previously  deposited  metal  and  is 
directly  below  the  new  bead.  The  ratio  of  remelted  metal 
to  the  new  metal  deposited  for  various  weld  processes  is 
shown  in  Table  III-2.  A  typical  macrograph  from  which  this 
type  of  data  was  extracted  is  pictured  in  Figure  III-4. 

/*-• 

\  «• 
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LATERAL  BASE  PLATE  MELT 
(INCH) 


WELD  PROCESS 


HEAT  INPUT 
(KJ/IN) 


1/32 

SMAW 

55 

1/32 

SAW 

55 

1/32 

GMAW  (PULSE) 

110 

TABLE  III-l 

LATERAL  BASE  PLATE  MELT 
FOR  VARIOUS  WELD  PROCESSES 


RATIO  OF  REMELTED  METAL 

TO  MEW  WELD  METAL 
DEPOSITED 

WELD  PROCESS 

HEAT  INPUT* 
(KJ/IN) 

0.6 

GMAW  (SPRAY) 

55 

1.0 

SMAW 

55 

1.0 

GMAW  (PULSE) 

110 

1.5 

SAW 

55 

TABLE  I II- 2 

REMELTED  METAL  TO  NEWLY  DEPOSITED 
WELD  METAL  RATIOS 


*Heat  input  of  55  KJ/in  is  the  current  maximum  allowed  by 
welding  specifications.  The  110  KJ/in  is  considered  high 
heat  input  and  was  done  for  experimental  purposes  only. 
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o  Analytical  Results 


In  1979  the  SSN704  weld  joint  was  re-analyzed  using  the 
B0S0R5  computer  code  considering  elastic-plastic  material 
behavior,  1/32  inch  lateral  base  plate  melt,  various 
remelt  depths,  and  a  stress  relief  force  system  to  simu¬ 
late  the  effects  of  re-melt.  For  treating  the  re-melt 
phenomenon  in  the  analytical  model;  where  the  weld  metal 
is  assumed  to  be  deposited  in  layers,  equations  1,  2,  and 
3  were  derived  to  calculate  membrane  and  bending  forces 
acting  at  the  mid-plane  of  the  weld  segment  equivalent  to 
the  stresses  existing  in  the  metal  which  is  assumed  to  be 
remelted.  The  negative  of  this  force  system  is  applied  to 
the  analytical  model  to  simulate  the  effect  of  relieving 
the  stresses  in  the  re-melted  metal.  See  Figure  III-2. 
Equation  1  AXIAL  LOAD 

fcR 

N  =  T~  (°Z1M  +  °Z2M} 


Equation  2 


MOMENT  LOAD 


M  =  a„ 


Z1  fcR 
t  (  — -  +  — 
■R'2  6  ' 


t  ('1  +  -*) 
k  '  2  3 


See  Figure  III-2  for  definition  of  terms. 
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The  results  of  the  B0S0R5  elastic-plastic  analyses  for 
1/2,  1,  and  2  layer  assumptions  of  remelt  are  shown  in 
Figure  III-3  together  with  the  measured  distortion  curve. 
Distortions  for  various  remelt  depths  were  predicted  in 
order  to  determine  the  effect  of  remelt  depth  on  distor¬ 
tion.  The  1/2  layer  remelt  condition  is  the  closest 
approximation  (about  10%  error)  of  the  SSN704  data  which 
is  to  be  expected  since  that  joint  was  welded  with  the 
GMAW  process. 

o  Summary  of  Results 

1.  For  the  SSN704  circumferential  hull  butt,  the  most 
recent  analytical  method  (BOSOR5  elastic-plastic 
including  1/32"  lateral  base  plate  melt,  single 
bevel,  AT  =  1300°F  and  1/2  layer  remelt)  yielded  a 
distortion  prediction  within  10%  of  the  measured 
distortion.  This  was  not  as  accurate  as  the  dis¬ 
tortion  predicted  by  the  1978  BOSOR5  elastic-plastic 
analysis  without  remelt  which  was  within  5%  of  the 
measured  distortion.  However,  both  predictions  are 
considered  within  engineering  accuracy. 


2.  Predicted  distortion  when  remelt  is  considered  via  the 
applied  stress  relief  force  system  increases  as  the 
depth  of  remelt  increases.  This  is  just  the  opposite 
of  what  one  would  expect.  Greater  remelt  should  cause 
more  stress  relief  and  less  distortion. 

3.  The  most  recent  analytical  method  which  includes  stress 
relief  is  fairly  sensitive  for  changes  in  remelt  pen¬ 
etration  into  the  previously  deposited  weld  metal.  As 
the  remelt  depth  was  increased  from  1/2  to  1  to  2 
layers  the  predicted  distortion  increased  from  .117" 

to  .149"  to  .182"  or  about  55%  over  the  range  con¬ 
sidered. 

4.  Based  on  weld  test  sample  evaluation,  a  base  plate 
melt  of  1/32"  on  each  side  of  the  weld  joint  is  the 
correct  value  for  use  in  the  analytical  model. 
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BEVEL  ANGLE 
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FIGURE  ffi-  I 


I* ROOT  8AP'«  MEASURED  ROOT  GAPt  BASE  PLATE  MELT 


FIRST  LAVER 


WELD  t  REMELTED 
METAL  ASSUMED 
TO  COOL  THRU 
AT«  -  1300  °F 


RE  MELTED 
LAYER 


LAYER  .  t 


. - -'  -FORCE  SYSTEM  DEVELOPED 

k(r,1M  _  BY  COOLDOWN  OF  WELD 
METAL  FROM  BOSORS 


NEGATIVE  FORCE  SYSTEM 
APPLIED  TO  RELIEVE  STRESSES 
IN  RE- MELTED  METAL 


BASE  PLATE  MELT 


LAYER 


/''*-~-____AVERAGE  MEASURED  BEVEL  ANGLE 

WHERE  ;  ~  _ _ _ _ _ _ 

Cl = (ROOT  GAPV  01)2 

C2  =  [TANGENT  (  BEVEL  ANGLE/ 2 )]  *  2 
H  =  RADIAL  LINE  LOAD 
M  =  LINE  MOMENT  LOAD 
N  =  AXIAL  LINE  LOAD 
R  s  MID  SURFACE  RADIUS  OF  SHELL 
tR=  THICKNESS  OF  WELD  LAYER 

Z,  =  RADIAL  DISTANCE  FROM  WELD  t  TO  INNER  SURFACE  OF  REMELTED  LAYER 
Z2  =•  RADIAL  DISTANCE  FROM  WELD  L  TO  OUTER  SURFACE  OF  REMELTED  LAYER 

crz ic  'CIRCUMFERENTIAL  STRESS  AT  DISTANCE  Z, 

<Tj|m*  MERIDIONAL  STRESS  AT  DISTANCE  Z, 

0"i2C»CIRCUMFERENTlAL  STRESS  AT  DISTANCE  Z2 
**2*' MERIDIONAL  STRESS  AT  DISTANCE  ?2 


FIGURE  JR  -  2 


WELD  BEAD  MACROGRAPH 
OF  SAW  WELD  TEST 
PLATE  WT  595 
(t3  *  MAGNIFICATION) 


(FULL  SCALE) 


FIG^E  IE -4 


IV  CONCLUSIONS  AND  RECOMMENDATIONS 


o  Heat  Transfer  Analysis  of  HY-130  Cylinder  Girth  Weld 
The  comparison  of  the  predicted  and  experimental  temperatures 

for  the  first  weld  pass  show  a  large  deviation,  up  to  200°F, 

for  a  surface  node  located  1.0  inches  away  from  the  weld  pool. 

The  amount  of  deviation  decreases  for  nodes  further  away  from 

the  pool.  The  probable  reason  for  this  discrepancy  is  the 

presently  uncertain  value  of  the  specific  heat  of  the  weld 

metal  in  the  liquid  state,  as  explained  in  Section  II-A-3. 

This  hypothesis  may  be  proven  by  a  two-step  approach,  as 

follows . 

1)  Perform  additional  computations  utilizing  trial 
values  of  the  liquid  specific  heat  until  the 
predicted  and  experimental  temperatures  agree. 

A  resultant  specific  heat  of  the  order  of  five 

or  less  times  the  solid  state  specific  heat  would 
indicate  the  likelihood  of  the  hypothesis.  An 
extremely  high  specific  heat  would  negate  the 
hypothesis . 

2)  Presuming  that  (1)  above  indicates  a  likelihood, 
determine  the  specific  heat  and  other  thermo¬ 
dynamic  properties  of  the  liquid  state  experimentally. 
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o  Shrinkage  Force  Method  of  Analysis  of  HY-130  Cylinder 


Girth  Weld  via  NONSAP 

The  results  of  the  subject  analysis  look  quite  good  by 
inspection.  Distortion  during  heat-up  and  cooldown  cycles 
follows  a  logical  pattern.  However,  when  comparing  to  the 
experimental  results,  it  was  found  that  for  the  first  layer 
of  weld  the  predicted  distortion  (.0054"  inward)  was  in  the 
opposite  direction  of  both  the  measured  and  calculated 
values  reported  in  Reference  3.  Interestingly,  an  examin¬ 
ation  of  the  residual  strain  magnitudes  measured  in  the 
test  seem  to  support  a  distorted  shape  as  indicated  in  Fig¬ 
ure  II-C-12  of  this  report.  Also,  the  transient  temperature 
distribution  measured  in  the  test  shows  that  the  HY-130 
cylinder  was  not  allowed  to  cool  back  down  to  the  200°F  pre¬ 
heat  temperature  where  predicted  and  experimental  distortions 
could  be  compared.  The  only  correlation  with  test  data  is 
the  initial  outward  distortion  shown  in  Figure  II-C-11  when 
the  weld  boundary  nodes  are  in  the  first  heat-up  cycle. 

A  trip  was  made  to  MIT  to  obtain  any  additional  test  data 
that  might  help  resolve  the  above  discrepancies.  It  was 
learned  that  no  distortion  or  temperature  test  data  was  avail¬ 
able  for  welding  after  the  first  bead.  An  examination  of 
strips  cut  out  of  the  HY-130  cylinder  in  way  of  the  girth 
weld  showed  that: 
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1.  The  first  weld  bead  did  not  penetrate  through 
the  bevel  land. 

2.  The  mating  bevel  lands  were  touching  (after 
welding) . 

Both  of  these  tend  to  support  the  Electric  Boat  Division 
boundary  conditions  and  model  at  the  plane  of  symmetry  (weld 
centerline).  However,  if  the  bevel  lands  were  not  clamped 
tight  together  (i.e.,  slightly  separated)  during  first  pass 
welding,  an  outward  deflection  as  measured  could  be  expected. 
Another  test  similar  to  that  performed  at  MIT  is  required  to 
resolve  the  discrepancies.  However,  to  obtain  maximum 
utility  from  the  experimental  data,  this  next  test  must  be 
tightly  controlled  and  well  documented.  Transient  and  final 
(cooled  down)  measurements  of  temperature,  distortion  and 
stress  are  required  for  every  weld  bead  and/or  weld  layer. 
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o  Submarine  Hull  Cylinder  Girth  Weld  Analysis 

The  results  obtained  when  applying  the  Shrinkage  Force  Method 
to  the  SSN704  hull  butt  raised  two  important  questions: 

1.  Why  did  the  distortion  predictions  incorporating 
remelt  stress  relief  loadings  increase  with 
remelt  depth? 

2.  Which  variation  of  the  Shrinkage  Force  Method  in 
the  axisymmetric  case  will  accurately,  cost- 
effectively,  and  consistently  predict  weld  distortions. 

o  Recommendations 

The  following  tasks  are  recommended  for  future  investigation: 

1.  In  the  heat  transfer  problem  perform  additional 

computations  varying  the  liquid  weld  metal  specific 
heat  until  predicted  and  experimental  temperatures 
agree  in  order  to  prove  or  disprove  the  hypothesis 
that  the  probable  reason  for  the  discrepancy 
between  the  predicted  and  measured  temperatures  is 
the  uncertain  value  of  liquid  specific  heat. 
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2.  Experimentally  determine  the  specific  heat  and 
possibly  other  key  thermodynamic  properties  of 
liquid  weld  metal. 

3.  Conduct  a  sensitivity  study  of  the  temperature- 
dependent  characteristics  of  the  mechanical 
(structural)  properties  of  the  weld  filler 
material.  This  will  identify  the  structural 
parameter  (or  parameters)  that  are  key  proper¬ 
ties;  and,  therefore  identify  those  areas 
where  high-temperature  property  testing  should 
be  concentrated. 

4.  Future  effort  should  also  be  devoted  to  the  develop¬ 
ment  or  acquisition  of  a  computer  program  for 

the  prediction  of  the  three-dimensional  temper¬ 
ature  distribution  in  the  weld  pool  and  base 
metal  including  the  effect  of  the  movement  of 
the  weld  rod.  It  is  potentially  possible  that 
transient  temperature  distribution  can  be 
obtained  from  the  equivalent  steady  state 
solution  of  the  problem  cast  so  that  the  thermal 
field  stays  in  a  fixed  position  relative  to 
the  weld  rod. 
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5.  Perform  a  tightly  controlled  and  well  documented 
test  similar  to  that  reported  in  Reference  5 
measuring  distortions,  temperatures  and  induced 
strains  for  each  weld  bead  or  layer. 

6.  Obtain  additional  hull  distortion  measurements 
in  way  of  submarine  circumf erential  butt  weld  to 
expand  and  increase  the  confidence  level  in 

the  experimental  data  base.  This  type  of  data 
is  invaluable  in  analytical  correlation  efforts. 

7.  Apply  the  various  formulations  of  the  Shrinkage 
Force  Method  to  the  submarine  hull  butts  measured 
in  6.  above  in  order  to  determine  which  formulation 
is  the  most  accurate,  cost-effective  and  consistent 
approach . 
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